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The dynamics of the tilted axisymmetric Bianchi IX cosmological models are explored allowing 
energy flux in the source fluid. The Einstein equations and the continuity equation are presented 
treating the equation of state w and the tilt angle of the fluid A as time dependent functions, but 
when analysing the phase space w and A are considered free parameters and the shear, the vorticity 
and the curvature of the spacetime span a three-dimensional phase space that contains seven fixed 
points. One of them is an attractor that inflates the universe anisotropically, thus providing a 
counter example to the cosmic no-hair conjecture. Also, examples of a realistic though fine-tuned 
cosmologies are presented wherein the rotation can grow significant towards the present epoch but 
the shear stays within the observational bounds. The examples suggest that the model used here can 
explain the parity-violating anomalies of the cosmic microwave background. The result significantly 
differs from an earlier study, where a non-axisymmetric Bianchi IX type with a tilted perfect dust 
source was found to induce too much shear for observationally significant vorticity. 


I. INTRODUCTION 

Cosmological observations provide compelling evidence that our Universe at large scales is well described by the 
homogeneous, isotropic and spatially flat standard model called the ACDM [1]. The common explanation for this 
featurelessness is that inflation, while also at the same time giving birth to small initial seeds of structure, wipes out 
all inhomogeneities from the exponentially expanding background. Furthermore, even the initial structure generated 
by inflation seems to be of a vanilla nature, as the data remains consistent with primordial spectrum from quantum 
fluctuations of a single field with perfectly gaussian and isotropic statistical properties [2]. 

However, not all the observational data comply with this picture, in particular the anomalies in the cosmic microwave 
background (CMB) [SHI]. If of cosmological origin, these anisotropic statistical features in the temperature fluctuations 
of the CMB could be regarded as a hint of physics beyond the simplest ACDM parameterisation of cosmology. As the 
anomalies are most significant at the largest CMB angles [5] , thus (roughly) corresponding to the scale of dark energy, 
it is natural to consider (slightly) imperfect fields in the present universe that could both accelerate the background 
expansion and generate the observed deviations from isotropy [31], see also e.g. HMD- On the other hand, the 
anisotropies could have been generated (or retained) by non-standard inflationary dynamics due to for example vector 
fields [T51E] or quadratic curvature corrections to gravity [23123| . 

However, such an anisotropic expansion during either the early (inflation) or the late (dark energy) accelerated 
stages of the universe, has a limited promise of actually explaining the CMB anomalies, because the latter are mainly 
parity-violating but shear does not distinguish handedness. It is this oddness that appears the most unexpected and 
difficult to explain with the existing models. Proposed late-time origin for a parity violation include spontaneous 
isotropy breaking [22], imperfect dark energy field [2S| and effects due to our peculiar velocity with respect to the 
CMB |26j . The primordial spectrum itself could contain dipole and other odd contributions only in the context 
of noncommutative theories ED, but parity-violating primordial modulations of the perturbations and thus of the 
resulting CMB temperature field might be generated by nontrivial multifield dynamics during inflation |28j or reheating 
[29j ■ or by breaking translation invariance with isocurvature perturbations |30j or domain walls [3I1E2I. 

The question we shall pursue here is whether our universe could be realistically described by a spacetime with slight 
anisotropies that distinguish orientation in way that would explain the parity-violating CMB anomalies. A natural 
framework for this is provided by the Bianchi class B models that describe homogeneous but anisotropic spacetimes 
that allow for rotation [221 [21]. Such are sufficiently general to incorporate parity-violating cosmological features 
(unlike Bianchi class A models that allow only for even asymmetries) but due to remaining symmetries (we consider 
locally rotationally symmetric Bianchi models) simple enough to be tractable. Indeed, it is has been shown that the 
CMB anomalies can be very well fitted by a Bianchi class B type Vllh template [21135], but only the anomalies: 
the cosmological parameters of the best-fit Bianchi template model are in strong disagreement with the rest of the 
cosmological data [251 [26] . This suggests it may be worthwhile to explore different or more elaborate Bianchi type B 
cosmologies. 

We are thus lead to consider generalisations of the above cited anisotropic inflation and imperfect dark energy 
models that have been previously discussed, as far as we know, only in the context of Bianchi class A spacetimes. 

A rather generic property of perfect-fluid cosmological spacetimes is that they tend to isotropise, which is also the 
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reason for the difficulty of reconciling the Bianchi Vllh rotating universe with realistic cosmology. This property is 
reflected in a cosmic no-hair conjecture suggesting that anisotropies will not survive in inflating spacetimes |37L I38j . 
Behind the conjecture is a theorem stating that initially expanding homogeneous cosmological models, which satisfy 
Einstein’s equations with positive cosmological constant and dominant and strong energy conditions, can evolve only 
so that they either collapse or approach the de Sitter solution at the same time making the universe isotropic. The only 
exception is Bianchi IX type universe, where the curvature can come with the positive sign and could thus in principle 
cancel the driving effect of shear and energy density to the expansion EZIEH]- As mentioned above however, several 
examples have been found where anisotropies can be supported in the presence of imperfect fluids, and refinements 
of the cosmic no-hair theorem have been formulated |3H| . 

Motivated by these considerations, we will specialise here to cosmologies in the Bianchi class B type IX models. In 
order to possibly generate vorticity, we include a ’’tilted” energy source, i.e. a fluid whose rest frame does not coincide 
with the cosmic rest frame, and furthermore it turns out that this requires that the fluid has an energy flux that is 
fixed by consistency conditions. When analysing the phase space, for simplicity, we will assume that the fluid can be 
characterised by a constant equation of state and that the tilt angle remains constant during the cosmic evolution. 
This defines our two-parameter phenomenological model of tilted dark energy in a universe allowing rotation. From 
this simple starting point emerge already interesting possibilities such as anisotropic inflation. Since Misner’s classic 
mixmaster paper |40j , studies of Bianchi IX cosmologies in the literature have been undertaken in the contexts of loop 
quantum cosmology [HI 02] , chaotic dynamics |43l 04] and viscous fluids 06] , but not to our knowledge in the 
presence of tilt and heat flux. 

The structure of the paper is as follows. We first specify the model in sectionj^and derive its governing differential 
equations, then analyse its phase space in section [Till in particular locate the fixed points there and investigate their 
stability properties. In section jlV] we then consider the possible cosmological applications of the model, focusing on 
the new anisotropic features we could generate without running into conflict with observational upper bounds on 
them. The paper is concluded by section jV] and a few details are confined to appendices |X]|C] 


II. THE MODEL 

We shall work in the framework of the axisymmetric Bianchi IX cosmological model, where the line element is [dOj 

ds^ = -d^ + e2“(‘)+2/3(b [(^1)2 + (^2)2] ^ g2a(t)-4/3(i)(^3)2 ^ 

with the one-forms w* given as 

oj^ = sin ijjd9 — sin 9 costpd(j), ( 2 ) 

= cos iIj d9 + sin 9 sin tpdcj), (3) 

oj^ = cos 9d(f) -I- dtp ■ (4) 

There are 5 non-zero components of the Einstein tensor using the line element Q, however, El = and 
El = cos(9)(Ejl — El), thus only three of them are independent. 

We consider the energy-momentum tensor to be composed from perfect fluid and energy flux, 

Tb = [pit) + pit)]u°'Ub + p{t)gl + q°-Ub + u°-qb , (5) 

where p{t) is the matter density, pit) is the pressure, and the energy flux vector qa is orthogonal to Ua- We take the 
four-velocity to be 


= (cosh[A(t)],0,0,sinh[A(t)]e2^(*)-“(*)) , 


( 6 ) 


which satishes the normalisation condition u‘^Ua = — 1. The ’’tilt” vanishes in the limit A —>■ 0. The fluid we assume 
to obey, 

pit) = wit)pit). (7) 

Using the non-diagonal components of the Einstein equations^ = Tip) and the orthogonality condition (gait® = 
0), the components of the energy flux vector are found: 

<Zo = -g3tanh[A(t)]e2«‘)-«d), ^^=0, g2 = g3Cosd, 93 = -(p(t) + p(t)) ^ (8) 

cosn[.2A(^t jj 


1 


We set SnG = 1. 
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Hence, after the dynamics of a, (3, p, p, A, and 9 are known, the dynamics of the energy flnx is describable. Therefore, 
we can study the dynamics of the system without paying attention to the energy flux and consequently we neglect 
it from now on. However, the presence of the energy flux in the energy-momentum tensor is crucial, because in the 
absence of the energy flux the equations would be consistent only if w = — 1 or A = 0. In a sense, the role of the heat 
flux is to allow w and A vary freely. The energy conditions for the energy-momentum tensor ([^ are not affected by 
the tilt and heat flux as discussed in the Appendix]^ Moreover, by defining an auxiliary function 


Wit) ^ I 


^1 - wit) + 


wit) +1 A 

cosh[2A(t)]/ ’ 


(9) 


the energy-momentum tensor ([^ can now be recast as 


T = 


-Wit)pit) 0 
0 wit)pit) 
0 0 

0 0 


0 

0 

wit)pit) 

\Wit) — l]p(t) cos(0) 


0 

0 

0 

[wit) + Wit) - l]pit) 


( 10 ) 


We can now determine the suitable kinematical scalars and spatial curvature of the metric, in terms of which to 
present the Einstein field equations. The covariant definitions of the kinematical quantities in the 3-1-1 formalism 
read as follows [34l 02] (where the covariant time and space derivatives are defined, respectively, as / = u“/a and 
f,a = i5l+U^Ua)f,b): 


Wab — ^[a;b] ; ^ab — ^(a;b) , 0 — U U — U . (H) 

These are, respectively, the vorticity, the traceless shear, the expansion and the acceleration tensors associated to the 
fluid and subject to the conditions uj(^ab) = o'lab] = 0: WabU^ = cTabU^ = 0i = 0 and = 0. From these quantities 
we can then construct the kinematical scalars, that for the line element Q and the four-velocity (§ turn out to as 

a; = (cCaf,a;“V2)'/' = ^sinh[A(t)]e-(“(*)+4'3(‘)), 

a = ((Ta&cr“'’/2)^^^ = A/3cosh[A(t)]/3'(t), (12) 

ii = iuaU°')^^'^ = sinh[A(t)][a'(t) — 2/3'it )], 

0 = 3cosh[A(t)]a'(t), 


where w is the vorticity, a is the shear, 0 is the expansion, and u is the acceleration scalar, respectively. The spatial 
curvature scalar (i.e. the Ricci scalar of the 3-dimensional constant-t hypersurfaces) for the line element 0 is 






-2(a{t)+im) 


(13) 


Using Eqs. ® , Q and (10) and the above definitions, the Einstein equations = T// assume the form of a constraint 
equation and two dynamical first-order differential equations: 


-Wit) pit) 
wit)pit) 
[wit) + Wit) - 1] pit) 


, ^jt)^ 

2 3cosh^[A(t)] cosh^[A(t)] 

a;(t)^ Qit)^ 0(t)CT(t) o'(^)^ 20'(t) cr'(<) 

sinh^[A(t)] 3cosh^[A(t)] -\/3cosh^[A(t)] cosh^[A(t)] 3cosh[A(t)] -s/S cosh[A(t)] 

^Rjt) , 2wit)^ ejt)^ 20(t)a(t) 

2 sinh^[A(t)] 3cosh^[A(t)] -\/3cosh^[A(t)] 

ajt)^ 2Q'it) 2a'it) 

cosh^[A(t)] 3cosh[A(t)] -y3cosh[A(t)] 


The conservation of stress energy WaT°'^ = 0 results in the differential equation 


0 = ^ (p(t)sech[A(t)] (eit)i3wit) + 4Wit) - 1) - 2y3E(t)(lU(t) - 1)) + 3 (p(t)lU'(t) + lU(t)p'(t))) 


(15) 
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However, we still need to determine the evolution of and uj in order to form an aut ono mous system. This can be 
achieved by differentiating both sides of the definitions of w and given in Eqs. (12)-(13), and then reapplying Eqs. 
pl-pl, giving 


= -^R{t) 


20(t) 


2tT(t) 


3cosh[A(t)] v^cosh[A(t)] 

0(0 , . cr(0 


4\/3a;(t)^cr(t) 
sinh[A(t)]2 cosh[A(t)] ’ 


3cosh[A(t)] v^cosh[A(t)] 


(16) 

(17) 


This completes the system of dynamical equations in terms of the variables in the 1+3 formalism [MlllZ]- 

To apply the techniques of dynamical system analysis in cosmology^ [33], it is convenient to rewrite the system 
in terms of dimensionless, expansion-normalised variables. Eor this purpose we find it most suitable to employ the 
Hubble rate that in the presence of tilt is defined by 


H{t) = 


0(0 


3cosh[A(t)] 


The dynamical variables we then define as follows 


m = 


r(0 


\/3 cosh.[\{t)]H (t) 


m = - 


QHit) 


2 ’ 


v(0 = 


lit) 


sinh[A(t)]i7(t) ’ 


n{t) = 


p(0 

3i7(t)2 


(18) 


(19) 


These are identified, respectively, as the dimensionless expansion-normalised shear, spatial curvature and vorticity 
associated to the spacetime geometry, and finally the energy density assoaciated to the fluid source. By assuming the 
universe to be monotonically expanding and treating a as the time variable, we have for any function fit) that 


f(0 /'(O 

H{t) ct'{t) 


= /'(«)■ 


( 20 ) 


This choice of time variable enables us to eliminate the Hubble rate from the resulting system of equations. Indeed, 
after some simple manipulations, Eqs. (14), (15), (16), and 0 can be written using the dimensionless quantities 
([T^ as 


S'(a) = -1 + 2K{a) + E{a)^ + (e - 3)E(a) + V{a)'^ + Vt{a ), 


( 21 ) 


H'(a) = H(a) (^2e + 2S(a) - 4 + 


l-3u;-2E(a)A H(a) 


W(a) 


W{a) 


W'{a) 


( 22 ) 


V'(a) = V(a)(e-4S(a)-l) , (23) 

Kfa) = -2 {K{a)[-€ + 1 + S(a)] + S(a)V(a)2} , (24) 

where we have defined e = —H'{t)/H{t)^. The four dynamical equations are subject to the constraints 

e = 1 [3w - 1] n{a) - K{a) + 2 + S2(a), (25) 

l = W{a)nia) + K{a) + T,‘^{a). (26) 


The first constraint means simply that the e in the four dynamical equations should be considered as a short-hand 
notation for the function of the four variables (the physical meaning of e is explained shortly), but the second equation, 
representing the Friedmann constraint, reduces the dimensionality of the system to three. 

Having now set up the dynamical system in a convenient form, we proceed to analyse its phase space. 


3 For more recent applications to anisotropic cosmologies, see for example the references [Jl [S] 1151117112111441 in the introduction. 
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III. THE PHASE SPACE ANALYSIS 


In this section, w and A are treated as variables for simplicity and we investigate the phase space of the three- 
dimensional autonomous system^ obtained by eliminating the curvature using the generalised Friedmann equation 
(26), employing the analysis methods described in many text books, e.g. [531135]. Even though the curvature is no 
longer among the set of differential equations, it can be solved using the generalised Friedmann equation. We will do 
this in the subsequent to track the curvature. 

We also follow e at each fixed point. That is of interest because this so called slow-roll parameter gives us directly 
information about the overall expansion rate of the universe. Namely, the geometrical mean, a, of the scale factor is 
a = e“, and hence. 




a 


H{t)^ (a')2 a • (a')2 

and thus e < 1 indicates acceleration and e > I deceleration. From Eqs. 


+ 1 1 


and (261 we furthermore see that 


= I 
= 3 

= 2 - 


1 -I- sech(2A) 


in the curvature filled universe , 
in the shear filled universe , 

in the dust filled universe, 

in the radiation filled universe, 
in the dark energy filled universe . 


At each fixed point in the phase space, e is a constant. 

We note that the tilt appears in the dynamical equations only in the combination W = {—w -I- 1 -b co™h[ 2 A) )/^- 
In the limit of vanishing tilt, W always reduces to unity. In the limit of strongly tilted cosmology, W reduces to 
W —>■ 1/2 for dust and W —> 1/3 for radiation. If ic = —1, the dynamical equations are explicitly independent on A, 
but the equation for the vorticity, (23), still holds. Thus, to possibly rotate the universe we need only a tilted fluid. 


The implicit dependence on A is that equation (23) does not exist iff A = 0. 


After these preliminary considerations, we are ready to list our results for the fixed points and their properties. 


A. The fixed points of the system 

We found seven solutions to the set of equations V' = ft' = T,' = 0. In the following, we will present these fixed 
points and study their stability in the parameter space (A,w), where —oo < A < oo and —1 < re < 1. We use an 
asterisk on quantities evaluated at each respective fixed point. 


® We could of course write down the three dynamical equations | |2l|22|23|24[ l with the constraints | |25|26| l inserted explicitly but that 
would hardly give more insight to the properties of the system. 

^ The governing equations depend on A only via cosh(2A). Due to the propeties of the hyperbolic cosine function, the studied range 
—oo < A < oo can be presented at the range 0 < A < oo. 
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W 



cosh(2 A) 



FIG. 1: The left panel shows the stability properties of the fixed point l |27| l: the orange area denotes saddle fixed points (repulsive in one 
direction) in the parameter space, the dotted curves are where the stability of the fixed points can not be determined using the linearisation 
method and in the gray zone V takes imaginary values. The vorticity and the curvature are zero on the curve between the orange and the 
gray areas. The right panel shows the corresponding values of e. Colors available online. 


1. The fixed points allowing rotation 

The two fixed points, 


3 {2 - 3w)w + 8{W - 1)W + 1 

^2y (1-5W)2 

3ii2-3w)w + 8{W-l)W + l) 

4(1-5W)2 ^ ’ 

3{w-6W + l) 

(1 - 5W)2 
3w + 2W-l 
low -2 

differing only by the sign of the vorticity V*, exist when the inequality 

0< (2-3w)w + 8(W-l)lV + l (28) 

is satisfied. It is convenient to note that condition ( [2^ implies that 0 < 51V — 1 must hold too. The area restrained 
by ( [^ is depicted in Figure where also the stability of the fixed points is indicated by colour. The left panel of 
Figure is drawn by evaluating the real part of each eigenvalue at each point separately. The eigenvalues are given 
in Appendix The dotted curves represent the parameter values for which at least one eigenvalue of the stability 
matrix is zero and the stability can not be determined using the linearisation method. Moreover, K* = V* = 0 on 
the dotted curves (for further details, see Appendix |B 1[ ). The right panel of Figure [^indicates that these fixed points 
can represent the overall expansion of a curvature, shear, dust, or radiation dominated epochs. 


V* = 
K* = 

n* = 

E* = 


2. The fixed point where only the vorticity does not exist 

The fixed point where all the effective energy sources contribute to the expansion (but the vorticity disappears) is 


V* 

K* 

fl* 


* 


0 

3 / (2 - 3w)w A 

4 ((1-2IV)2 ^ J 

3w 

“(1-21V)2 
3u; + 21V - 1 


(29) 


E^ 


41V- 2 
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cosh(2 X) 


cosh(2 A) 



FIG. 2: The left panel shows the stability properties of the fixed point ( |29| l. In the orange areas the solution is a saddle point (repulsive 
in one directions in the lighter colored regions and in two directions in the darker colored regions), in the white area the fixed point is 
stable, and on the black curves the fixed point does not exist. Dotted curves (including u; = 0) denote where the linearization method 
can not recognise the stability properties. If A ^ oo, then 1 and w^, and converges to zero. The right panel shows the 

corresponding e-values. Colors available online. 


The stabilities and values of the slow-roll parameter e are depicted in the parameter space in Figure These fixed 
points can present any expansion rate as e can take the value of any real number. The curve where the fixed points 
are not defined is w = [cosh(2A) — 1]”^, furthermore e —> oo when approaching the curve from the negative w side 
and e —> —oo when approaching from the positive w side. This feature is absent without tilt. In the region where the 
fixed points are stable, i.e. when 

< w < 0, (30) 

we can have 1 < e < 3/2. For the eigenvalues and curves wf, see Appendix 

The full stability of the w = 0 solutions can not be solved using the simplest Lyaponov method. The linear method 
showed the fixed point to be stable on this curve in the shear and the vorticity directions. Using the center manifold 
theorem, we find the lowest order approximation in the dust direction to be 


B2 


n'{a) 


2sech^(2A) 
sech(2A) + 1 




(31) 


From this equation can be deduced that the fixed point is unstable in the dust direction: if the system is perturbed 
in the positive dust direction, it will return to the fixed point, but if the system is perturbed in the negative dust 
direction, it will continue in the negative direction. Consequently, when w = 0, the system has no stable fixed points. 


3. The fixed point allowing stable inflation with shear-hair 


The nontrivial fixed point where the vorticity and the curvature disappear is at 


V* 

K* 

n* 

T,* 


0 

0 

3(u>- l)(3w-4IU + l) 
W{-iw + 2W + IY 
2-2IF 


(32) 


—iw + 2W + 1 


This situation is depicted in Figure]^ The fixed point can represent any expansion rate as it exists for —oo < e < oo. 
The curve where the fixed points are not defined is ic = (1 + 2cosh(2A))/(—1 + 4cosh(2A)), moreover e —oo when 
approaching the curve from the negative w side and e —>■ —oo when approaching from the positive w side. This feature 
is absent without tilt. There exists a wide region of parameter space where the fixed point is stable, i.e. when 


w < w_ 


B3 


The eigenvalues and w_ are presented in Appendix 
w « — 1. Therefore, this fixed point gives a counter example to the cosmic no-hair conjecture. 


Furthermore, the right panel of Figure shows 


(33) 

e « 0 when 























cosh(2 A) 


cosh(2 A) 



FIG. 3: The left panel shows the stability properties of the fixed point | |32| |. In the white region the fixed point is an attractor, and a 
saddle point in the orange regions (repelling in one direction in the lighter and two directions in the darker shaded areas). On the black 
curves the fixed point is not defined, and on the dotted curves (including w = 0) a Lyaponov exponent has a vanishing real part. If A ^ oo, 
then 1, tCg 1/2, 1/5 and 0. The right panel shows the values of e. Colors available online. 


4 . The shear universe fixed points 

The system includes two fixed points where only the shear do not vanish: 

= 0 and E* = ±1. 

The eigenvalues of the negative dimensionless shear fixed point are 


w 

and the eigenvalues of the positive dimensionless shear fixed point are 

-iw + AW-l 


- 2 , 2 , 


W 


(34) 


(35) 


(36) 


Both fixed points are always saddle and have e = 3 for all w and A. 

5. The curvature and shear universe fixed points 

The system includes another fixed point independent of the tilt angle and the equation of state parameter: 

V*=C* = 0 and X*=3/4 and S* = 1/2. 

The eigenvalues of the fixed point are 


iw 

'W 


(37) 


(38) 


These are all negative if w S (0,1], thus describing a stable fixed point. The solution is not an attractor if the fluid 
dilutes slower than dust, w € [—1,0), since then relative contribution of the fluid energy density will dominate in the 
end. In the case where w = 0, the stability can not be determined using the linear method, but as in this case this 
fixed point coincides with (29l, we already know it is unstable. This solution expands with the rate given by e = 3/2, 
regardless of w and A. 


IV. COSMOLOGICAL APPLICABILITY 

In this section we consider the possible cosmological relevance of the models. First we will discuss the role of 
dynamical variables of the system (19), and then construct explicit examples of early and late cosmologies with 


nontrivial but potentially viable anisotropic dynamics. 
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A. On the anisotropies 

Here we briefly discuss various possible effects that seem to emerge in these anisotropic models. In particular, we 
report the observations that 1) the shear-free condition could be supported with the tilt 2) the rotation might grow 
in special cases 3) strong tilt can have a freezing effect on the evolution and 4) we obtain anisotropically accelerating 
solutions. The last two of these we illustrate with some numerical solutions obtained by integrating our full equations 
motions. The numerical solutions are in complete agreement with the analytic considerations of the previous section. 


1. Shear 


Cosmological observations constrain the shear to be very small in the post-recombination universe. In particular, 
the CMB is compatible with at most |S| < 10“'^ at the decoupling [5]. This appears to be problematic if we want to 
have otherwise sizable anisotropies, since according to Eq. (HI), curvature as well as vorticity act as sources for the 
shear. Thus, generically in viable cosmologies we expect to have only small effects from each of these terms, since it 
is otherwise difficult to keep the shear under control. We will present in the following some examples with nontrivial 
but small cosmological effects from anisotropies in the presence of shear. We have found both stable and unstable 
fixed points with nonvanishing shear, of which (321 can accommodate accelerating expansion and can thus describe 
anisotropic inflation. 


With special fine-tuning however, one might eliminate the shear from the anisotropic fixed points (27) and (291 (the 
shear is zero in the fixed point ( |32[ ) iff A = 0 or w = — 1). The fine-tuning required is the specific relation between the 
fluid parameters cosh2A = —(1 -I- l/u>)/2 that singles out a zero-measure family of models in the w < 0 half of the 
parameter space. We will not investigate these special cases further here, but mention that it might be interesting 
to consider them as a possible realisation of the so called shear-free condition in cosmology, see Refs. [Hiiniisn] on 
theoretical and phenomenological aspects of shear-free cosmologies. 


2. Rotation 


The vorticity determines an infinitesimal rotation in an infinitesimal proper time interval m- Its role differs from 
the shear and curvature anisotropies, as the vorticity decouples from the expansion. By this we mean that it does 
not enter the constraints ( 26|25 ) and thus does not act as an effective energy source for the expansion, but rather 
affects the dynamics of the overall scale factor indirectly via its dynamical impact on the evolution of the curvature 
and the shear. Furthermore, the redshift of a distant object is not directly dependent on V either, because it does not 
affect on the separation between two successive spatial hyper surfaces m- The fixed point where the vorticity is not 
zero, (271, is always saddle and can represent the shear, the curvature, the radiation or the matter dominated epochs 
according to possible e values. 

From the evolution equation (23) we see that the condition for the vorticity to grow is e > 1 -|- 4E. Thus, in a 
shear-free case we cannot increase the rotation in an accelerating universe. Negative shear could however catalyse the 
growth of vorticity. 


3. Tilt 


For strongly tilted models, cosh(2A) 3> 1, the time passing along the trajectories of the system may slow down. Let 
us illustrate this ’’freezing” effect with numerical solutions. Choosing a positive-density matter with the equation of 
state re = 0, we experimentally find that the fixed point (37) appears to attract a wide range of solutions. The rate 


at which the solutions converge to this attractor however can depend considerably upon the tilt angle A of the matter 
fluid. This is seen in Fig. where trajectories approaching to the fixed point ( |37[ ) are drawn with appreciable 
tilting angle. We set the tilt to A = 3 and evolve models with various initial conditions, all the curves appear to stop 
evolving after 5 e-folds, where V « 0 and E « 0.5 (the left panel). However, allowing the integration to continue 


much longer, say 200,000 e-folds longer, we notice that all the curves do still approach the fixed point (37) along the 
line V = 0 and S = 0.5 (the right panel). In contrast, with say A = 0.1, the system is close to the fixed point (37) 


already in tens of e-folds. One of the dynamical effects of the tilt is thus obviously that it slows down or ’’freezes” the 
evolution (at least in some corners of the phase space near the fixed points). 
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FIG. 4: Trajectories in a dust-filled universe with A = 3. In the left panel, the system has evolved for for 1000 e-folds starting from 
different initial values, but none of the trajectories seemingly has evolved after 5 e-folds. In the right panel, the same trajectories have 
had additional 200000 e-folds time to evolve, in which time all of them have nearly reached the fixed point Pit . Hence, after 5 e-folds the 
trajectories do evolve slowly along the same line parallel to the fl-axis, where V = 0 and S = 0.5. See SectionjlV A 3|for further discussion. 


4- Inflation 


The fixed point (291 in the system can represent a universe in a stable exponentially expanding phase (e « 0). The 


shear is non-zero at the fixed point, unless A = 0 or w = — 1. If w is slightly larger than negative unity, the universe is 
approximately in an inflationary phase regardless of the tilt. However, a non-zero tilt makes the expansion faster and 
induces shear. Hence, this fixed point represents a counter example to the cosmic no-hair conjecture. As expected, 
this requires violationg of the strong energy condition (namely, w < —1/3, see Appendix]^. The amount of shear is 
connected to the rate of expansion in this point, and is thus not arbitrary if we want to have inflation. In the limit 
of large tilt, these quantities are given solely by the equation of state of the fluid: 


1 -b w 

2 — 4w 


3 (1 + w) 
2 — Aw 


when A —>■ oo . 


(39) 


These limiting values are reached with good accuracy with A of order few. Remarkably, we can have both the shear 
and the expansion rate within the desired ranges. Namely, the near-de Sitter inflation with e ^ 0.05 that is required 
to produce the observed tilt of the scalar fluctuation spectrum [2] is perfectly compatible with a shear at the level of 
few percents, which is the maximal amount allowed by constraints on the anisotropy |35j (in fact, with large enough 
tilt simply S* = e/3). 

Numerical examples of cosmologies converging to the anisotropic inflationary attractor are shown in Fig. with 
the models parameters set as w = —0.99 and A = 4. Trajectories are shown with four qualitatively different initial 
conditions, including an isotropic case corresponding to standard inflationary expansion. We see that, from all the 
initial conditions, the universe ends up in the inflationary fixed point that efficiently brakes the rotation but retains 
a finite shear supported by the tilt. 


B. Bianchi IX cosmology with late-time anisotropies 


We will construct here an explicit example of full cosmological evolution. The purpose is to find a concrete scenario 
wherein the cosmological parameters fit the observations and the anisotropies are not too large to be immediately 
ruled out by the CMB bounds, but wherein we could still generate a small amount of rotation (that could have 
relevance to the large-angle CMB anomalies). We consider a scenario where dark energy is tilted. The equations to 
integrate are Eqs. (2l]|26), but in order to describe full realistic cosmic history, we need include also dust and radiation 
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FIG. 5: Four solutions of the system, where w = —0.99 and A = 4, are drawn in the vicinity of the stable fixed point | |29| l. There the 
system is in exponentially accelerating state, e ^ 0.005, where V* = 0, S* Ri 1.7 X 10“^ and 12* Ri 1.005, hence presenting a counterexample 
to the cosmic no-hair conjecture. The solid curve starts from where an isotropic inflation would occur, but it also evolves to the anisotropic 
fixed point. 


in the play and supplement the corresponding terms in these equations. At low redshifts we then expect that our 
tilted dark energy component will begin to dominate the energy budget (over dust and radiation) and in addition to 
accelerating the universe, introduce slight anisotropic features. The aim of the numerical study is to investigate the 
detailed dynamics within the latter transition period as well as at the time of decoupling. 

In this set-up then dust and radiation are comoving with the coordinates and the dark energy fluid is tilted with 
respect to them. To distinguish dust, radiation, and dark energy from each other, we can still treat the variables (19) 
as effective ones and merely decompose the energy-momentum tensor ([^ as 


= {P2 + P2)n°"n^ + + (P3 + P3)n°'n^ + + (Pi + Pi)uiu\ + + QiA + ulq\, 

where = ( 1 , 0 , 0 , 0 ) is the four velocity comoving with the coordinates and 

< = (cosh(Ai),0,0,sinh(Ai)e2'^(‘)-“(‘)^ 


(40) 


(41) 


is the four velocity of the tilted fluid and Ai is the exact tilt of the tilted fluid with respect to the comoving coordinates. 
The details of the decomposition are given in Appendix [C| 

By using the energy-momentum tensor ( |40[ ), the relevant field equations become 

e = i [3w3 - 1] fl3(a) + ^ [3w2 - 1] f^ 2 (a) + ^ [3wi - 1] Oi(a) - K{a) + 2 + , (42) 


1 = n^ia) + fl 2 (a) + ITifli(a) -b K{a) + E^ia ), 


(43) 


when it comes to the two constraint equations, and the four dynamical equations of motion generalise to 
= —1 + 2K((y.) H- — 3)S(gi) H- V(o;)^ H- fli(Qi) H- ^^ 3 ( 0 ;) 5 


(44) 
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n[ia) 


^ , 1 — 3wi — 2E(a) 

f 2e + 2S(a) — 4 H- — - 

fi 2 (a) + ll 2 (a)[ 3 (l + W 2 ) — 2 e] + 113 (a) + ll 3 (a)[ 3 (l + W 3 ) — 2e] 

iVi 


(45) 


V'(a) = (e-4S(a)-l)V(a), 


(46) 


K'(a) = -2 {K(a)[-€ + 1 + E(a)] + E(a)V(a)2} , (47) 

where IFi = (—'*i'i + l +“ Kpi/{3H'^) and =Pi/pi [i = {1,2,3}). Like in Sectionpl the conservation 
of energy-momentum tensor yields only one equation, even though the decomposition of results two new functions. 
To obtain a closed system, we are required to impose two new equations. The geometrical mean of the expansion is 
e“, hence imposing pi oc 5 - 3 ( 1 +“^)“ g^d in terms of the dimensionless variables 


r2'(a) = r2i(a)[2e - 3(1 -I- w^)] 

Using this relation for 112 and II 3 , Eq. (451 simplifies to 

1 - 3w - 2E(a) 


ll'i(a) = lli(a) (^2e-g2S(a)-4-g- 


Wi 


(48) 


(49) 


After choosing appropriate values for the parameters Wi and Ai, Eqs. (42), (43), (44), (46), (47), (48) and (49) 
constitute a full autonomous system which we use to describe the evolution of the universe from the radiation 
dominated epoch to the present time. Eventhough earlier we discussed about a set up where dark energy is tilted 
with respect to dust and radiation, the aforementioned equations are derived without assuming which fluid is tilted. 
Hence in the subsequent, we give special attention to tilted dark energy scenarios, but employ other possibilities too. 


1. Numerical solution 


We studied the above specified system of equations (42), (43), (|44|, (46), (47), (48) and (49) numerically in order to 
quantify the amount of fine-tuning required for the system to describe viable late time cosmo ogy, including the CMB 
anomalies. In practice we mean by viability now that with anisotropies kept at the maximum of some percent level, 
the system should evolve through radiation and matter dominated eras as usual, and then near the present epoch 
approach a dark energy dominated solution so that IIa ~ 0.7 and fid « 0.3 and fir ^ 10“® today. The time from the 
decoupling to the present time then corresponds to In(llOO) « 7 e-folds (with shear at the maximum percent level, 
the expansion is direction-independent to a good approximation). From Eq. it appears that the tilt shouldn’t be 
too strong in order to avoid generating too much shear in the expansion. We thus impose A < 1/100, and furthermore 
choose the fluids to be dark energy, dust and radiation. 

According to our numerical analysis, the solutions satisfying the above described conditions all originate from the 
vicinity of the negative or positive fixed point (34). However, not all the initial conditions close to the fixed points 

From the vicinity 
but can lead to different late 


(34) lead to the desired type of solutions, but the initial conditions need to be carefully chosen, 
of either of the fixed points (34), the solutions enter the radiation dominated era alike, 
time cosmology. Nevertheless, constraining K and V from above appears to be sufficient to keep the anisotropies 
negligible, independently of the sign of E, which fluid is tilted or the amount of tilt (as long as A < 1/100). These 
results indicate that inflation is not required to vaporise shear for viable late time cosmology to occur. Let us consider 
examples where some shear is left after inflation. 

Example A: The tilted fluid is dark energy and A = 1/100. The initial conditions are given in Table [l] and the 
solution is drawn in Fig. (§, showing the dynamical variables as functions of the number of e-folds (a). In addition, 
the values of the dynamical quantities at the decoupling and at the present time are tabulated in Table |l] We observe 
the emergence of three epochs: the radiation, the dust and the dark energy dominated expansions, the latter finally 
converging to the attractor (29) as expected. Also the shear, the vorticity and the curvature are drawn in the figure. 


though they are almost indistinguishable. They reach the peak of their dynamical significance around the matter - 
dark energy equality, where we roughly have E « —0.0007, K « —0.004, and V « 0.05. Eventually these will decay 
with the accelerated expansion, since also shear given by the asymptotic future solution (29) vanishes when ic = — 1. 


Let us study study this example further. According to our numerical simulations, at the beginning of the radiation 
epoch only conditions \K\ < 10“^° and |V| < 10“® are to be satisfied to ensure small enough anisotropies during late 
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time cosmology. These conditions are remarkably alike to those stipulated for the standard ACDM model referred 
as the flatness problem. Breaking the condition for K would make either |E| too large at the decoupling or both |E| 
and \K\ too large during the transition from dust dominated to dark energy dominated eras. Also the vorticity is 
constrained from above, since too much rotation would boost the E and K during the transition period. In addition to 
these initial conditions, of course, we need to choose the relative amounts of the three types of energy densities, fli’s, 
suitably to have proper radiation, dust and dark energy epochs. The initial condition for E is then to be determinate 
using Eq. (43). Further numerical analysis indicates that, at the accuracy presented in Table |T] and Figure]^ the 
results of Example A are the not dependent on the size of A (assuming A < 1/100) and hold for both signs of V and 
irrespectively which fluid is tilted. 

Example B: Let us consider the same scenario as in Example A, but change the sign of the initial shear. The 
initial conditions are given in Table |l] and the solution is drawn in Fig. <§• In addition, the values of the dynamical 
quantities at the decoupling and at the present time are tabulated in Table |Tj The system appears identical to that in 
Example A apart from E and V. From the figure can be seen that the ratio |V/E| takes the higher value at a ~ —0.5 
in the case of negative initial shear and at the same a the shears take opposite signs. From Table 1 we see that 
these conclusions hold at the decoupling and present time too. Further numerical analysis indicates that the same 
constraints as in Example A, namely \K\ < 10“^° and |V| < 10“^, apply here for the same reasons. Moreover, at 
the accuracy presented in Table |l] and Figure the system in this example does not alter by changing the fluid that 
is tilted, the size of A( < 1/100) or the sign of V. 

Example C: Let us study Examples A and B using different values for initial vorticity. We impose the initial 
values 


(0.01 X Vab, 0.1 X Vab, Vab, 10 x Vab), 


where Vab denotes the initial value for V used in Examples A and B, i.e. Vab = 1 x 10“^, and represent the results 
in Table [n] and in Figure]^ Furthermore, the evolution of E and V corresponding to the initial value X x Vab is 
denoted by Ejc and Vx in Figure respectively. For example, the curves Eq.qi and Vq.oi describe the evolution of E 
and V corresponding the initial vorticity 0.01 x Vab- The effects of changing the initial V appear alike in both cases. 
Curves Eg.oi and Eq.i are indistinguishable for both cases in Figure]^ indicating initial V < O.IVab has negligible 
effect on the evolution of E. Increasing initial vorticity further induces positive contribution to shear and vorticity in 
both cases around the time of dust - dark energy equality. Qualitatively the behavior is similar but quantitatively 
weaker at the decoupling and at the present time (see Table [ll|. Increasing the initial vorticity V < Vab by an order 
of magnitude induces an order of magnitude growth for V but negligible growth for E at the decoupling. The results 
are independent on the amount of tilt, assuming A < 1/100, and which fluid is tilted. 

The peculiar feature that all the results exhibit independence on the size of A < 1/100 and which fluid is tilted 
imply that the evolution of the anisotropies is decoupled from the fluid for small enough tilt.® In the corner of the 
phase space studied above, the anisotropies reach their peak values during matter - dark energy transition. This 
is no surprise, because after that point they can no longer grow as they approach the asymptotic future solution 
(29) cutting all hair of anisotropy when w = —1, hence this location is at least a local maximum of anisotropies. 


Because also shear and vorticity appear to be decoupled for a small enough tilt and vorticity at the beginning of the 
radiation era, the results obtained here differ significantly form another Bianchi type IX realisation |51j . where even 
small vorticity was found to induce considerable shear. In Ref. [51] the authors found the maximum value for the 
present time vorticity to be ~ 3.9 x 10“^® for non-axisymmetric Bianchi IX metric sourced by tilted perfect fluid 
including only dust. Around 10 magnitude difference in the results suggests that different realisations of the Bianchi 
IX space-time can yield remarkably different results. In the light of Ref. |51] and our numerical examination, the 
model used here shows potential explaining the anomalies in the CMB and late time expansion simultaneously. 


V. CONCLUSIONS 

We investigated the dynamics of the tilted axisymmetric Bianchi IX model with fluid allowing energy flux. The 
cosmological application we had in mind was the generalisation of anisotropic inflation and anisotropic dark energy 
models to generate parity-violating imperfect features at cosmological scales. The latter could perhaps then be 
manifest in the large-angle anomalies of the CMB. In the framework of the Bianchi IX model, accommodating a tilt 


® To support this conclusion, we explored the range 1/100 < A to a lesser extent and found the anisotropies begin to grow during the 
domination of the tilted fluid. The only exception was tilted dark energy, because w = —1 makes A disappear from the dynamical 
equations. 
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TABLE I: The initial values for the dynamical variables of the models in Examples A and B are given, the exact value of shear 
is to be calculated from the generalised Friedmann equation (431, but only its approximative value is presented in this table. 
The systems are solved using these initial data and the values of the dynamical variables at the decoupling time and at the 
present time are reported with one digit accuracy. 



Da 

Dd 

Q.r 

E 

V 

K 

Example A - positive initial shear 

Initially 

1 X 10"^® 

5 X 10"® 

0.99 

« 7.07107 X 10"® 

1 X 10“® 

-1 X 10“^° 

At the Decoupling ~ 

1 X 10"® 

0.7 

0.3 

9 X 10“® 

2 X 10"® 

-6 X10"® 

Present Time ~ 

0.7 

0.3 

1 X 10“'‘ 

-6 X lO”'* 

4 X 10"® 

-3 X 10“® 

Example B - negative initial shear 

Initially 

1 X 10"^® 

5 X 10"® 

0.99 

« -7.16378 X 10"® 

1 X 10“® 

-1 X 10“^° 

At the Decoupling ~ 

1 X 10"® 

0.7 

0.3 

-8 X10"® 

3 X 10"® 

-6 X10"® 

Present Time ~ 

0.7 

0.3 

1 X 

1 X 10“® 

7 X 10"® 

-3 X 10"® 


TABLE II: The values of V and S at the decoupling and at the present time corresponding to different initial values for 
vorticity are presented. Here, 0.01 x Vab, 0.1 x Vab, Vab, and 10 x Vab represent different initial vorticities (see Example 3) 
and on the left side the values correspond to positive initial shear (see Example A) and on the right side the values correspond 
to negative initial shear (see Example B). 



0.01 X Vab 

0.1 X Vab 

Vab 

10 X Vab 

0.01 X Vab 

0.1 X Vab 

Vab 

10 X Vab 


At the decoupling 

At the decoupling 

V 

3 X 10"® 

3 X 10“'‘ 

3 X 10"® 

3 X 10“® 

2 X 10"® 

2 X 10-^ 

2 X 10"® 

2 X 10"® 

E 

-9 X 10"® 

-9 X 10“® 

-9 X 10"® 

3 X lO-^ 

8 X 10"® 

8 X 10“® 

9 X 10"® 

2 X lO"'* 


At the present time 

At the present time 

V 

7 X lO"'^ 

7 X 10“® 

7 X 10"® 

3 X 10“® 

4 X lO""^ 

4 X 10“® 

4 X 10"® 

3 X 10"® 

E 

-2 X 10"® 

-2 X 10“® 

-1 X 10"® 

4 X 10“® 

-1 X 10"® 

-1 X 10“® 

-6 X 10“'‘ 

3 X 10"® 


of the fluid required us to introduce a heat flux for consistency. In general, the Einstein equations and the continuity 
equation are given treating the equation of state w and the tilt angle of the fluid A as time dependent functions, 
enabling us to consider realistic multi-fluid scenarios, but for simplicity the phase space was studied by reducing w 
and A to parameters. 

We performed a dynamical system analysis, identifying seven fixed points in the cosmological phase space and 
determining their stability for given model parameters w and A, which we considered in the ranges — 1 < ic < 1 and 
—00 < A < 00 . Only one of the fixed points allows for rotation, but always describes non-accelerating expansion. The 
fixed points generically include a non-negligible amount of shear, which can then be accompanied by the curvature 
with or without the tilted fluid energy density. Of particular interest for possible applications is an inflating attractor 
in which the amount of shear is given by the tilt of the fluid. The vorticity in general has a different status in the 
system compared with the other dynamical variables, as the vorticity does not contribute directly to the expansion 
rate. 

The anisotropically accelerating fixed point can be relevant to both early inflation and dark energy models. We 
did some preliminary numerical investigations of both types of scenarios, verifying that indeed realistic cosmological 
evolutions can be constructed. The late time solutions show potential explaining both the parity preserving and 
violating anomalies of the CMB, which is unexpected in the light of earlier studies of the subject [ST]. Our numerical 
studies indicate that the only fine-tuning the late time solutions require is to set the curvature very close to zero after 
inflation. The fixed point (29) retains a constant amount of shear (when w ~ — 1 and A yf 0) with the expansion, 
even if it is accelerating, thus providing a counter example to the cosmic no-hair conjecture. Though many such 
counter examples have been presented (see the Introduction |I] for an incomplete list of references), such have not been, 
previously to our knowledge, deviced by exploiting tilt and energy flux. 

It could be very interesting to study further the observational implications of such possible imperfections in the 
source fluid. Determining the precise signatures from relaxing the standard assumptions q°‘ = Q and A = 0 could 
allow us to falsify the proposed origin of cosmic anomalies as a tilted dark energy (or inflaton). However, before 
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FIG. 6: Bianchi IX cosmological models, Example A on the left side and Example B on the right side, describin g th e evol ut ion o f th e 
universe from radiation dominated epoch to the dark energy dominated epoch. The system is solved using Eqs. | |42[ |, l |43[ l, 

( |49| and l |48[ l and the initial conditions in Table 1 . The horizontal axis is the number of e-folds from the present time. All the 
evolving quantities are plotted in the upper panels, whereas the behavior of V, S, and K is magnified close to the horizontal axis in the 
lower panels. 



FIG. 7: The evolution of S and V in two Bianchi IX cosmological models is integrated using several initial values for vorticity, where 
the initial shear is positive (left plot) and negative (right plot). Curves denoted by Vi and Si represents the evolutions of the vorticity 
and the shear, respectively, and the subscript refers to the initial value of the vorticity (see Example 3). The curves Sq.oi and Sq.i are 
indistinguishable from each other in both panels. 


such proposal could be promoted from the present phenomenological exploration into a more convincing alternative 
cosmology, we needed to develop also the theoretical underpinnings of the model. It would be desirable to describe the 
fluid by a lagrangian field theory, instead of adding it as a parameterised energy component by hand. In particular, 
in the present context the precise form of the heat flux was fixed by mathematical consistency rather than physical 
motivations. To describe the dynamics of fluid more properly, we would guess one probably needs to allow in general 
time-dependent w(a) and A(q!), and the dynamics for these should be determined from the lagrangian principle for 
the fundamental degrees of freedom for the fluid. 

Despite these challenges, we believe our results could provide an interesting starting point for further explorations 
of rotating universes filled by some kind of tilted dark energy. 
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Appendix A: Energy conditions 

The energy-momentum tensor ([^ does not violate the weak or dominant energy conditions. For some parameter 
values w and A, it does however, break the strong energy condition. The strong energy condition for ([^ reads, using 
the results of Ref. [52] for our case, 

p ■y/(w -I- 1)2 tanh^(2A) -I- tc -f 1^ > 0 p {w + l)^sech2(2A) -I- 2w'j > 0 . (Al) 

For w > 0, the above conditions are satished for any A, but the latter of the conditions is violated by any A if 
w < —1/3. Thus neither the tilting of the fluid nor the implied energy flux changes the usual bound for the violation 
of the strong energy condition. 


Appendix B: Eigenvalues 


1. The rotating fixed point (27| 


The eigenvalues £/ of the stability matrix of the fixed points (271 are 
3w + 2W-l 


£■/ = 3 ■ 


5VF- 1 


^2 =3 


VW{w -W) + VA/2 
^/W{bW -1) ’’ 


^3 =3 


vW(w - kF) - VA/2 
vW(5kF- 1) 


(Bl) 


where A = —2w -I- re — 3) -I- 8(u> -I- 14)kF^ -I- 40('u; — l)wW — 92W^ — 28W + 2. The dotted curves separating the 

different areas in Figure are 


^ cosh^(2A) — 2 ^1 ± a/ 3 — 2 cosh(2A)^ cosh(2A) -b 1 


w+ = 


cosh(2A)(cosh(2A) -b 4) — 2 


(B2) 


Curves w/ satisfy the equality in (281 and equations £2 = K* = V* = 0. 


2. The non-rotating matter-scaling fixed point (291 


The eigenvalues of the stability matrix of the fixed points (291 are 


^2 3 3w -b 21F - 1 

1 “ ~2 2kF- 1 ’ 


£7 = -v 


3 W{w-2W+ 1)-VB 2 _ 3 “ 2kF-b 1)-b 

t'O —- 7 


(1 - 2W)W 


(1 - 2W)W 


(B3) 


where B = W {—2Aw^ -b w‘^{W -b 16) -b w (28kF^ — 30kF -b 8) -b (1 — 2kF)^kF). The curves separating the different 
areas in Figure]^ are obtained from equations £i = Q (dotted) and l/Ei = 0 (solid), i = {1, 2, 3}. These curves are 




-1 


Wn = 0 , w\ = 


1 2 cosh(2A) -b 1 

^ 2 cosh(2A) -b 1 ’ cosh(2A) — 1 ’ cosh(2A) — 1 ’ 

1 — cosh(2A) -b cosh^(2A) ± cosh(2A)y^cosh^(2A) — 2cosh(2A) -b 5 
2cosh^(2A)-b 2cosh(2A) — 1 ’ 


(B4) 


(B5) 


and it is easy to verify that they approach to zero or unity as given in Figure The -b case of Eq. (B5) yields the 
lower bound for inequality (30). 
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3. The non-rotating flat matter-scaling flxed point (32| 


The eigenvalues Ei of the stability matrix of the fixed points (32) are 

3(ri;-l)(3w-41T + l) i + 2w + {I - 2Wf) + Qw + 24W'^ - 2AW + 3 

“ —?rF 777 — , in —) ^2 — -7777—7— , -, 03 — - 


2W{-3w + 2W + 1) 


W{-3w + 2W + 1) 


-6wW + 41T2 -h 21T 


(B 6 ) 


The curves separating the different areas in Figure|^are obtained from equations £i = 0 (dotted) and l/£i = 0 (solid), 
i = {1, 2,3}. These curves are 


wf = 1 , W 2 = 


3 cosh(2A) + 2 3 


1 - 


5cosh(2A) - 2 ’ 


cosh(2A) + cosh^(2A) ± cosh(2A)'ycosh^(2A) — 2 cosh(2A) + 5 
2cosh2(2A) + 2cosh(2A)-l 


, (B7) 


_3 2 — 2cosh(2A) + cosh^(2A) ± 2cosh(2A)\/3 — 2cosh(2A) 3 2cosh(2A) + l 

^ cosh^(2A) + 4cosh(2A) — 2 ’ ^ 4cosh(2A) —1 

At the limit A —?► 00 , the behaviour of curves w^, W 2 , and wf is given in Figure]^ 


Appendix C: Decomposition of the energy-momentum tensor 


fluid, whose density consists 
The effective or the total equation of state parameter w 


In the set of Eqs. (21)-(26), the fluid variable ft can be considered as an effective ’’total' 
of a combination of separate components as = Eifli 
describes how the density of the total fluid is related to its pressure and the effective tilt illustrates how much the 
total fluid is tilted with respect to the normals of the surfaces of homogeneity. In this sense, the quantities E, K, V, 
and are also to be considered as effective. Therefore, we can decompose fl and q°‘ into the ’’original” components, 
but still treat the other quantities in their effective forms. 

Consider two four-velocities, uj = (cosh(Ai), 0, 0, sinh(Ai)/Y^ 533 ) and n'’ = (1,0, 0,0). The average (or effective) 
four-velocity, is the linear combination au^ + bn‘^, where a and b are constants so that u'^Uc = — 1. By introducing 

a shifted angle A implicitly via 


<.= 45hA, and i, = cosl,(A)-A(^cosh(A.), 


sinh(Ai) 


sinh(Ai) 


the new four-velocity is nothing but rt'” = (cosh(A), 0, 0, sinh(A)/^/^). Hence, the tilt angle A can be interpreted as 
an average (or effective) tilt of the combination of the fluids with respect to the reference frame. 

The relations between the effective and exact fluid densities, state parameters and tilt angles can now be obtained 
by equating tensors (H and (|40[). However, a simpler expression is achieved by expressing (|40[) using IFi, 


/ 


T = 




0 


2 = 2 




2=1 


'^Wi{t)pi{t) 


0 

0 

cos( 6 ')(VFi(t) - l)pi(t) 


\ 


2=1 


0 ki(o + Wi{t) - i]pi{t)+y^wi{t)p^{t) 


2 = 2 


(Cl) 


and equating (10) and (Cl I, yielding 


pV) = P^^t) + ft(() + p,(t), .»«) = + + W{t) = + ,C2) 


Pl{t) + P2{t) + P3{t) 


Pl{t) + P2{t) + P3{t) 


Substituting these expressions to Eqs. (21) - (26) yields Eqs. (42) - (47) if Wi and Ai are constants. 
































